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Abstract 

The magnetic phase diagrams, and other physical characteristics, of the hole- 
doped Lof2,- x Sr x CuO^ and electron-doped Nda- x Ce x CuOi high-temperature 
superconductors are profoundly different. Given that it is envisaged that 
the simplest Hamiltonians describing these systems are the same, viz. the 
t — t' — J model, this is surprising. Here we relate these physical differ- 
ences to their ground states' single-hole quasiparticles, the spin distortions 
they produce, and the spatial distribution of carriers for the multiply-doped 
systems. As is well known, the low doping limit of the hole-doped mate- 
rial corresponds to k = (§, f ) quasiparticles, states that generate so-called 
Shraiman-Siggia long-ranged dipolar spin distortions via backflow. These 
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quasiparticles have been proposed to lead to an incommensurate spiral phase, 
an unusual scaling of the magnetic susceptibility, as well as the scaling of the 
correlation length defined by £ _1 (a;,T) = £ _1 (z,0) + £ _1 (0,T), all consis- 
tent with experiment. We suggest that for the electron-doped materials the 
single-hole ground state corresponds to k = (tt, 0) quasiparticles; we show 
that the spin distortions generated by such carriers are short-ranged. Then, 
we demonstrate the effect of this single-carrier difference in many-carrier 
ground states via exact diagonalization results by evaluating S(q) for up to 
4 carriers in small clusters. Consistent with experiment, for the hole-doped 
materials short-ranged incommensurate spin orderings are induced, whereas 
for the electron-doped system only commensurate spin correlations are found. 
Further, we propose that there is an important difference between the spa- 
tial distributions of mobile carriers for these two systems: for the hole-doped 
material the quasiparticles tend to stay far apart from one another, whereas 
for the electron-doped material we find tendencies consistent with the clus- 
tering of carriers, and possibly of low-energy fluctuations into an electronic 
phase separated state. Phase separation in this material is consistent with 
the mid-gap states found by recent ARPES studies. Lastly, we propose the 
extrapolation of an approach based on the t — t' — J model to the hole-doped 
123 system. 
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I. INTRODUCTION: 



The Cu02 plane based high-temperature superconductors have anomalous normal state 
properties. One part of the normal state puzzle involves the spin orderings, and for the 
majority of this paper we focus on the differing magnetic phase diagrams for the hole- 
doped Bednorz-Miiller La 2 ^ x Sr x Cu0 4: compounds |T[ in comparison to the electron-doped 
Nd,2- x Ce x Cu04 materials ||. In particular, for both of these systems the x = phase 
possesses three-dimensional antiferromagnetic long-ranged order [EM. However, one must 
heavily dope the Nd compound to destroy this order, say ~ 12 %, while only a small doping 
of ~ 2 % is required to kill long-ranged antiferromagnetic order in the La cuprate material. 
A schematic of the contrasting phase diagrams is given in Fig. 1. 

One important aspect of the magnetic orderings found in these systems involves the 
kinds of spin correlations that these systems exhibit when doped to levels greater than 2 %. 
As mentioned above, for the electron-doped material antiferromagnetic order persists until 
~ 12% doping levels are reached. However, when the hole-doped system has say 7 % holes 
in a CuOi plane, short-ranged incommensurate magnetic ordering is found. This appears 
experimentally in the dynamic magnetic response, and has been found via inelastic neutron 
scattering for a number of intermediate doping levels 0-0. Consequently, an alternate way 
in which one can phrase the question concerning the differing magnetic phase diagrams is: 
why does the electron-doped system maintain commensurate antiferromagnetic order at 
doping levels that induce incommensurate correlations in the hole-doped system? 

Several theories have addressed the phase diagrams of these doped antiferromagnetic 
insulators. When the La-based system is doped, the carriers predominantly occupy sites, 
in particular as 0~ states, whereas for the iVd-based material the carriers are believed to 
be associated with the addition of an electron to a Cu 3d 9 ion, thus forming Cu + . Then, a 
common explanation of these differing magnetic phase diagrams follows from the assumption 
of the complete localization of the carriers, and the ensuing spin distortions generated by 
such static defects. Namely: (i) The dramatic reduction of the Neel temperature T N (x) 
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in the hole-doped system results from frustrating ferromagnetic bonds being embedded in 
an antiferromagnetic background, the so-called Aharony model ||. Such frustrated bonds 
generate a long-ranged spin distortion with dipolar symmetry, and numerical studies pfl 



have shown that a version of this perturbation renormalized by quantum fluctuations p, 10 



indeed alters the spin-spin correlation length in a fashion consistent with experiment |TT 



(ii) Fixed Cu + sites in a background antiferromagnetic lattice act like static vacancies, and 



subsequently diminish the spin-wave stiffness in a manner consistent with experiment 



|i"3fl . Monte Carlo studies of this quantum dilution problem fT^Jl^l do not suffer from the 



minus sign problem, and thus are capable of accurately characterizing the spin correlation 
length. 

Unfortunately, these cannot be complete explanations of the differing magnetic phase 
diagrams. Near the antiferromagnetic phase boundary (temperature vs. doping) the carriers 
are mobile in both of these systems |T5|JT6| , and thus one must come to an understanding of 
the differing reductions of the spin orderings for mobile not just localized carriers. This is 
one of the focuses of this paper |T7|] . 

Here we put forward a proposal explaining the electron-hole asymmetry of the magnetic 
phase diagrams shown in Fig. 1 for mobile carriers. Our work builds on our earlier study 
|18| that showed that for the t — J model and the hole-doped system, a tendency towards 



short-ranged incommensurate order vs. hole doping arises. The results of are consistent 
with experiment 0-0] , and in part is further verification of earlier exact diagonalization work 
of Moreo et al. as well as Monte Carlo work | 20| , pT |. However, the important point of 



T8| was that this incommensurability only occurred for the single-hole system having a 
ground state momentum of k = ±(|,±|) (since k = (f,f) quasiparticles generate long- 
ranged spin distortions similar to those produced by a ferromagnetic bond, the success of the 
above-mentioned Aharony model bodes well for any theory predicting that these carriers 
exist in the hole-doped 214 systems). If the momentum of the single-hole ground state is 
k = (it, 0), something that may be accomplished through the use of the t — t' — J model (see 
below) with t'/t > 0, then the tendency towards incommensurability is eliminated and only 



commensurate ordering remains. Clearly, this is similar to the above-mentioned behaviour 
of the hole and electron-doped materials, and in this paper we will make this comparison 
complete. 

Very shortly after ||18|| , Tohyama and Maekawa also studied the t — t' — J model, 



now for the electron and hole-doped systems. They suggested, consistent with a consider- 
able amount of electronic structure work p3|-|25||, as well as angle-resolved photoemission 
spectroscopy results for the electron-doped system , that the Nd system corresponded 
to hopping integrals t and t' with (approximately) just flipped signs in comparison to the 
same hopping integrals believed to be appropriate for the hole-doped materials (we elabo- 
rate on this relation in the next section). Then, by examining the energy to the first excited 
states, they demonstrated that the commensurate antiferromagnetic order was much more 
stable for the electron-doped system than for the hole-doped systems. As we show below, 



this is complimentary to our work in [Jl8| , since the single hole ground state momentum for 
their electron-doped t — t' — J Hamiltonian is actually k = (tc, 0)! We then strengthen the 
arguments of Tohyama and Maekawa by evaluating the magnetic structure factor S(q) for 
the many-carrier ground states for these two systems. Indeed, using the relevant material 
parameters, we find incommensurate ordering tendencies for the hole-doped material, while 
for the electron-doped material the magnetic ordering is always commensurate, in complete 
agreement with experiment. 

We shall also be concerned with other physical properties of the hole and electron-doped 
materials. For the La cuprates, some experiments can be explained beginning with the as- 



sumption that at low doping levels Shraiman-Siggia |27[ dipolar quasiparticles are present 



(implicit in these theories is the existence of hole pockets, or a "small" fermi surface, some- 
thing that so far has not been observed). One such magnetic problem involves the scaling of 
the correlation length with doping and temperature, viz. £ _1 (a;,T) = £ -1 (a:,0) + £ -1 (0,T), 
and this has been found to be reproduced by a model (based on the existence of these dipolar 
quasiparticles) proposed by one of us |2S| for both mobile and localized carriers, and is in 



close agreement with experiment [|TTJ. Further, a theory [[^| of the temperature and doping 



dependent magnetic susceptibility x( x i T) has been found to be consistent with experiment 
P0fl . Lastly, the same theory |H| also predicts the existence of a pseudogap, something that 
may be consistent with experiment ||32j| . 

We are proposing that these dipolar quasiparticles do not exist in the weakly-doped Nd 
cuprates, and thus, e.g., £(x,T), x( x ,T), and the optical properties of this system, may be 
quite different from those of the La cuprates. Unfortunately, the collection of experiments 
performed on the electron-doped materials is not as complete as the set performed for 
the hole-doped system, and so a full comparison is not possible. However, there is one 
striking difference between the normal state properties of these systems, and that is the 



resistivities at optimal doping. The La system shows a linear-T resistivity (T^|, whereas 
the Nd system shows a clear fermi-liquid-like T 2 dependence [Q. Another intriguing 
experimental result for the electron-doped material corresponds to the mid-gap states seen 



in the angle-integrated photoemission work of Anderson et al. ||26|| . It has been proposed 



34| that this is a signal of phase separation, and to this end we have examined the spatial 
distributions of many carriers, as well measured the interaction energies between carriers. 
Our results tend to support a phase separation hypothesis, but only for the electron-doped 
system. 

Our paper is organized as follows. In § |T| we discuss the formal aspects of the t — if — J 
model and display the electron-hole symmetries and asymmetries that can be identified when 
these materials are doped away from half filling. In § |T| we summarize previously published 
and our recent exact diagonalization results for the one-carrier ground states for these two 
systems, and show that these two systems have different one-carrier ground states. Then, 
the spin correlations found in S(q) as a function of doping are shown to be in agreement 
with the phase diagram of Fig. 1. In § [TV] we examine the spin distortions produced by each 
of these quasiparticles via our semi-classical analysis of this Hamiltonian for a single carrier; 
this will be shown to agree with the exact diagonalization results. Quantum fluctuations are 
shown to not affect these conclusions. Then we examine the spatial distribution of carriers 
pointing out the acute difference between the hole and electron-doped systems — from this 



we relate the tendency towards clustering of carriers of the electron-doped system to the 
mid-gap states seen in photoemission. Finally, in § [V] we discuss the relevance of our results 
to the transport properties of both kinds of systems, and as a prelude to possible future work 
we consider the extrapolation of this approach to other systems, and in particular consider 
the commensurability that persists in the hole-doped YBaCuO system. 



II. t — t' — J MODEL: 
A. Electron Hole Asymmetry induced by t': 

We begin our discussion of the t — t' — J model by first considering a simpler model 
involving only hopping terms to the nearest and next nearest neighbour sites along a one- 
dimensional chain — this model provides the most direct demonstration of the electron-hole 
asymmetry about half filling that if introduces into the problem. 

Consider a linear chain with a single orbital per site, and define the vacuum to be the 
state with all orbitals devoid of electrons. Introduce 

Hhop = ~ t ( C la C jo- + h.c)j - if ( C L C ka + h.c)j (2.1) 

<i,j>nn d <i,k> nnn d 

where < > nn denotes nearest neighbours, < > nm next nearest neighbours, and c ia is 
the electron annihilation operator for site % and spin a (note that these operators are not 
constrained to prohibit double occupancy — see § |11 U| ) . Then imagine that the system has 



an up spin at every site except one, at which no electron is present. The energy eigenvalues 
for this hole plus spin polarized background are given by 

E(k) = -2tcos(k) + 2t'cos(2k) (hole). (2.2) 

Now consider a system with up spins at every site, and that in addition to these electrons one 
down-spin electron is placed in some orbital. Then the energy eigenvalues of this composite 
spin-polarized plus mobile electron system are 

E(k) = -2tcos(k) - 2t'cos(2k) (electron). (2.3) 



From these results it is to be noted that while the near-neighbour hopping leads to 
identical energies for both the hole and electron-doped systems, the next-nearest-neighbour 
hopping has different signs. This reflects the breaking of electron-hole symmetry that t' 
induces, and thus its inclusion in our starting Hamiltonian is crucial. It is to be stressed 
that this asymmetry is found in any dimension, for any spin background, but only for doping 
levels close to half filling. 



B. Strong Coupling Hamiltonian: 



Equation (|2.1|) is not a good starting point for describing the physics of the C11O2 planes 
in either LaSrCuO or NdCeCuO, since it only describes free carriers. If instead one begins 
with an extended Hubbard model description |35| of the strong correlation effects known to 



be present in these systems, and then examines the systems at close to half filling, one can 
obtain the so-called t — if — J model ||36|| . To be specific, one considers 

h — — t J2 (4^ + h - c ) - E (44. + h. c ) + j £ ■ Sj - \. 

(2.4) 



<i,j>nn cr <i,k>„ nn a <i,j>nn 



where J is the Heisenberg superexchange constant coupling spins Si between nearest- 
neighbour sites. Also, in Eq. Q2.4|) constrained electron operators are used (constrained 
to disallow double occupancy at any site), viz. 

Cicr = Qct(1 c\_ (J Ci— cr ) . (2-5) 

We wish to stress that our use of the t — t' — J model is based on beginning with a 
three-band extended Hubbard model description []35] for both the LaSrCuO and NdCeCuO 
systems. While it is well known [p^-|39| that for the hole-doped systems one may approx- 
imately map the three-band Hubbard model onto the one-band t — t' — J model (N.B. - 
not necessarily the one-band Hubbard model), it may be argued that this is also the correct 
starting point for the electron-doped systems. For example, the dominant contribution to 
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the t! hopping integral for the electron doped system is not due to a direct Cu — Cu overlap, 
but rather follows from a third-order Cu — ► O — > O — ► Cm hopping |j4"0| . Thus, in what fol- 



lows we consider (i) the t and t' parameters that are to be fitted via quantum cluster studies 
(see below), and not due to a strong-coupling mapping of the one-band Hubbard model onto 
a Hibert space prohibiting double occupancy [^TJ, and (ii) the superexchange constant J to 



be found from either quantum cluster studies, or from comparison to experiments probing 
the magnetic structure and/or excitations. 

In the previous subsection we were able to treat hole and electron-doped systems with 
the same starting Hamiltonian. However, now it is clear that the above Hamiltonian, for a 
vacuum corresponding to no electrons at any site, cannot be used to describe both hole and 
electron doping of the half-filled state. For example, a vacancy can easily be added to a spin- 
polarized ferromagnetic state, and is completely mobile, while a down spin electron added 
to such a spin background is completely localized. Consequently, with this Hamiltonian and 
this vacuum, only hole-doped Cu0 2 planes may be investigated. 

To investigate the electron-doped system it is usual to introduce an electron-hole trans- 
formation, viz. 

Cia -> 4x (2.6) 

where is the annihilation operator of a hole at site i of spin a. (With this transformation 
it follows that one now considers a vacuum with no holes at any site.) However, this does not 
lead to the desired Hamiltonian as one finds that the new hopping term does not describe 
constrained (viz. two holes are never allowed on the same site) hopping. Instead, one must 
utilize a more general model than Eq. ( |2.4|) , one that includes so-called doublon hopping 



processes |[42|| , and from this starting Hamiltonian it is straightforward to derive the required 
transformation. The result of this transformation can then be shown to be equivalent to 

Cia -> 4t ( 2 - 7 ) 

where the also satisfy fermionic anti-commutation relations. Thus, one eventually finds 



that the appropriate t — t' — J model for the electron-doped system is the same as Eq. ( p.4| 
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but with the a^ replacing the q CT and with the minus signs in front of t and t' changed to 
be plus signs. 

Summarizing the above discussion, the t — if — J Hamiltonian will be used to model 
the physical systems that we are interested in. We consider this model to possess the bare 
minimum number of processes through which both the hole and electron-doped C11O2 planes 
can be represented. To be specific, (i) Eq. Q2.4D with a vacuum of zero electrons at every 
site of a square lattice will describe the hole-doped system, and (ii) for the electron-doped 
systems we use 

H = t (ataj* + h.c/j + t' E (SiffOfcff + h.c/j + J & ' ^ ~~ J n i n j) 

(2.8) 
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For this latter Hamiltonian we use a vacuum of zero holes at every site of a square lattice. 
Thus, note that, quite simply, in comparison to the hole-doped Hamiltonian the signs of 
the t and t' hopping terms are flipped and the a operators replace the c operators. Thus, 
from now on we simply use Eq. Q2.4| ) for both systems, specifying t = 1 for hole doping, and 
t — — 1 for electron doping . 



We now consider the numerical values of the hopping and superexchange parameters. A 
consequence of the holes residing on the O sites in the La material is that for a quantitative 
determination of the hopping parameters t and t', something very different from overlap 
integrals must be evaluated. We refer the interested reader to the summary of this problem 



given by Tohyama and Maekawa p2[ , and simply state the range of accepted values p2|-p4|. 
For the hole-doped LaSrCuO system, using Eq. (|2.4j) , one may scale all energies such that 
t = 1. Then one has J ~ .3 — > .4 and t' ~ —.2 — > —.4. For the electron-doped NdCeCuO 
system we again use Eq. ( |2.4| ) but now with t = —1, and J ~ .3 —>■ .4 and t' w +.2 — > +.4. 
Thus, it is interesting to note that even though a Zhang-Rice singlet is not exactly equivalent 
to a Cu 3d 8 ion, the material parameters for these two systems are such that the electron- 
hole mapping given in Eq. ( ]2.7| ) predicts reasonably well the numerical values that one must 
employ for one system given the parameters of the other. 
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III. NUMERICAL RESULTS: 



Stephan and Horsch jKJ have proposed that knowledge of the single-hole system does 



not lead to valuable information on the multiply-doped state. In [18] we showed that this 
is not always the case, and that, in particular, knowledge of the one-hole ground state can 
aid in predicting the presence/absence of incommensurate correlations. Thus, we begin a 
discussion of our numerical work with a summary of the single-hole ground state. 

There has already been a comprehensive exact diagonalization study of the ground state 
of one hole in the t — t' — J model for a 4 x 4 square cluster in Ref. H4J (since they use a 
bosonic representation, their sign of t' is flipped in comparison to ours). For the physically 
relevant ratio of J ft they find a k = ±(|-,±|-) ground state when the t' corresponding to 
hole-doped materials is used (unless | t' | becomes too large), and a k — ±(vr, 0), ±(0, n) 
ground state for the t' associated with the electron-doped materials. 

We have used a Lanczos routine and have also performed exact diagonalization evalua- 



tions |45| of the ground states of the t — t' — J model for J = .4, with t = ±1 for t' = ±.1, ±.2 
and ±.3. For the conventional 4x4 16-site square cluster our work agrees with that of Ref. 
44fl. We have also studied the a/8 x a/32 16-site cluster introduced in fl8fl, as well as some 



results generated using a V4~8 x a/32 24-site cluster P6] . For the hole-doped material, viz. 
when t' < 0, we always find a ground state of k = ±(f , ±f ), while for the electron-doped 
material, viz. when t' > 0, we always [§8| find a ground state of k = (±7r, 0), (0, ±7r). Our 
results for these two different clusters are entirely consistent with the phase diagram of Ref. 



44| found using the 4x4 cluster. 

Our numerical results also lead us to expect that the band structures for the hole and 
electron-doped materials will be quite different. We find that the band structure of the 
hole-doped material is similar to Fig. 4 of |18[ (which is a t' = curve), although with the 
increased ratio of | t'/t \ that we are using here the band along k = (0, 0) — > (tt, 0) becomes 
quite flat. This is similar to recent photoemission work of the valence band of undoped 



Sr 2 Cu02Cl 2 and provides some experimental support for a value of t' around —0.3 t. 
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We have listed the single carrier band structure energies for our y/18 x 24-site cluster 
in Table I for both the hole and electron-doped systems having used if /t — — 0.3. 

The single-hole ground states of the hole and electron-doped materials are different - 



so what? As mentioned above, the important fact contained in Ref. is that a knowledge 
of the single-hole ground state can aid in understanding whether or not any short-ranged 
incommensurate correlations are introduced into the ground state of the multiply-doped 
system. This followed from (i) if the single-hole ground state is not a k = ±(f , ±f ) state, no 
tendencies towards incommensurate correlations occurred, consistent with the spiral phase 
prediction of Shraiman and Siggia |50| , and (ii) the electron momentum distribution function 
for 2, 3 and 4 holes was found to be composed of the half-filled fermi surface with dimples 
at those momenta corresponding to the single-hole ground states, thus suggesting that some 
form of rigid band filling [18|] is in effect. Combining this with the differing single-hole ground 
states mentioned above, it is clear that we can now extrapolate these abstract results to the 
physical systems under consideration. 

Figures 2 and 3 show the magnetic structure factor for 0, 1, 2, 3, and 4 holes (electrons) 
for if = ±.3 evaluated with the v^8 x a/32 cluster introduced in [0 - similar results are 
obtained for other ±t' pairs. For the hole-doped material, Fig. 2 displays that the peak in 
the structure factor shifts from (ir,ir) for and 1 hole, to (^, ^f) for 2 holes, to (|, |) for 
3 and 4 holes. This behaviour is similar to the behaviour displayed in Fig. 5 of Ref. for 
if = 0. However, in Fig. 3 it is seen that no shifts of the maximum of S(q) away from (jr, tt) 
occurs for any doping concentration away from half filling for the electron-doped system. 
These results are entirely consistent with Fig. 1, and, in particular, are consistent with the 
lack of any tendency of the electron-doped materials to display the kind of incommensurate 
correlations that are present in the hole-doped system [|5T|. Clearly, these numerical results 
suggest that we have found a reasonable starting point from which one can hope to be able to 
describe the hole and electron-doped 214 systems. We now try to come to an understanding 
of these numerical results. 
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IV. ANALYSIS: 



As mentioned in the introduction, the existence of Shraiman-Siggia quasiparticles in 
weakly doped LaSrCuO is supported by the theories P, p8| , p9| , |3T| , |52| that assume their ex- 
istence, and are subsequently able to produce results consistent with experiment. To be 
specific, it seems that a scenario in which these quasiparticles are weakly interacting and 
tend to stay very far apart from one another (thus making it likely that some form of rigid 
band filling (e.g. with hole pockets around ±(f , ±f )) is applicable. 

Now we study the (71", 0) quasiparticles: we are proposing that these are the single-hole 
ground state constituents that exist in the weakly doped NdCeCuO system. We wish to 
characterize the spin distortions that each such quasiparticle produces, as well as come to 
an understanding of the spatial distribution for a state which is multiply doped. This will 
allow us to (i) understand why the (71, 0) and (|, |) quasiparticles lead to differing magnetic 
structure factors vs. doping, and (ii) examine some of the other physical features that are 
known to be specific to the electron-doped system. 



According to the semiclassical theory of Shraiman and Siggia p7|, for ground states with 
these wave vectors, differing distortions of the spin background occur. For the k = (| , | ) 
states, long-ranged spin distortions with dipolar symmetry are produced via backflow. For 
the k = (7r, 0) ground states, the hydrodynamic theory of Shraiman and Siggia [27| leads to 
the prediction that the spin distortion induced by one (ir, 0) quasiparticle is short ranged. 
These ideas are verified in our Figs. 4 and 5 where our results from evaluating the energies at 
this wave vector for the t — t' — J model with the electron-doped material parameters using 



a semiclassical variational wave function, such as that described in p7| , are displayed. To 
be specific, we utilize a product state of classical spins in an infinite lattice, and incorporate 
a broken antiferromagnetic symmetry in calculations of the hopping matrix elements for a 



single hole (see [^7] for the details of such a variational wave function). Then, to display the 
range of the spin distortions for this ground state we have calculated the minimum energy as 
a function of the number of spins away from the hole that are allowed to be distorted from 
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their undoped Neel configuration — these sites are shown in Fig. 4, where they are labeled 
according to the equivalent sites around a single carrier. In Fig. 5 we show the minimum 
energies found from the variational principle as a function of these site labels for both the 
hole-doped k = (|, ~) one hole ground state, and the k = (ir, 0) electron-doped one carrier 
ground state. The long-ranged spin distortion of the hole-doped ground state is clear, 
whereas the short-ranged distortion of the electron-doped ground state is made manifest 
by the insensitivity of the minimum energy when one allows more distant neighbours to be 
distorted from their Neel state. 

The above procedure may be repeated for all wave vectors for both hole and electron- 
doped Hamiltonians. We find that providing t' satisfies certain inequalities (which necessarily 
depend on the ratio of J/t) that the semiclassical band structure largely reproduces the 
band structures determined by exact diagonalization. Specifically, the ground states are 
reproduced when (i) t' > —.12 for the hole-doped material, and (ii) t' > +.12 for the 
electron-doped materials (these inequalities will be affected by quantum fluctuations — see 
below). This agreement with the quantum cluster studies provides further evidence that the 
renormalized classical description of the spin background (not including the charge carriers) 
of the weakly doped systems is valid at least at low temperatures H. 



This semiclassical work ignores quantum fluctuations. However, as shown by Reiter |53 
for an infinite lattice when quantum fluctuations are included in such variational wave func- 
tions one still finds the Shraiman-Siggia dipolar quasiparticles with the long-ranged dipolar 
spin distortion intact. We have found [|54j that within this lowest-order self-consistent ap- 



proximation |53| the (tt, 0) states still only involve short-ranged spin distortions, the same 
as the semiclassical prediction. Thus, it is not unexpected that our quantum cluster and 
semiclassical analyses agree. 

Now let us consider the effect of the short-ranged spin distortions that this kind of 
quasiparticle produces, viz. why do these quasiparticles not introduce any incommensurate 
correlations into the spin texture. Given the short-ranged nature of the spin distortion, 
it is clear that they could lead to a disturbance of the spin texture similar to that found 
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in the quantum dilution problem flllj, since static vacancies also produce a short-ranged 
disturbance of the spin texture. However, the quantum dilution problem begins with the 
assumption that the dilutants are randomly distributed throughout a plane. We now show 
that this is probably not the case for the mobile carriers in the electron-doped materials. 

In Ref. [^] we presented some numerics for the hole-hole correlation function for a pair 
of holes doped into a v8 x \/32 cluster described by the t — J model. We found that 
the holes tended to stay as far apart as possible. (It is interesting to note that this is 
entirely consistent with experimental work purporting to see charge-rich walls, separated by 
distances scaling like 1/x, arising from finite-size striped magnetic domains at low x [j55| . 



Thus, our work demonstrates a magnetic mechanism for such stripes, and is different than 
the Coulomb-interaction generated carrier-rich walls of Emery and Kivelson's [[34j phase 
separation ideas.) Since this numerical result is in direct contradiction to the results of 



Poilblanc ||56|| , who studied pairs of holes in a variety of square clusters, we have further 
scrutinized this behaviour. We have found that in (i) Monte Carlo simulations of the t — U 
model, for an average of 2 holes in an 8 x 8 lattice, the holes still tend to want to be as 
far apart as possible, and (ii) the same behaviour is found for the larger \^18 x \f%2 24-site 
cluster with 2 holes — these results will be published in a future comment [^7 . 



We now present our results for similar calculations for the electron-doped system. Fig- 
ure 6 shows the carrier-carrier correlation function, defined analogously to the hole-hole 



correlation function of Eq. (6.1) of ||18|| , by 



PJ\i-3\) = ttVE < (! - £«vO(l -E^) > 

where nj )(T is the hole number operator for electron-doped systems (electron number operator 
for hole-doped systems, as in Eq. (6.1) of [|18|]) for spin o at site i, Nt being the number 
of lattice sites, and N e is the number of equivalent sites a distance \i — j\ from site i (note 
that we do not include the angular dependence which is actually present for our non-square 
lattice, for simplicity only, since we have found that this does not affect our conclusions). 
We evaluated this function for two carriers doped into a \/8 x a/32 cluster described by the 
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t — t' — J model with t = — 1, J = .4 , and t' = +.3. Juxtaposed with this curve is the 
analogous hole-hole correlation function for t — 1, J = .4, and t' = — .3. The differences 
between the two different systems is striking: whereas pairs of holes tend to stay as far apart 
as possible, pairs of mobile electrons in the electron-doped system tend to cluster together. 

To further elucidate this behaviour we have calculated the 2, 3 and 4 carrier interaction 
energies (analogous to the binding energy) defined relative to independent carriers, (denoted 
by E21, E 3 i and En, where E n denotes the ground state for n carriers), via 

E 2 = E + 2(E 1 — E ) + E 2 i 
E 3 = E + 3(E 1 — E ) + E 3I 

E A = E + A{E 1 - Eq) + E AI (4.2) 

and our results for both hole and electron-doped systems are presented in Table II. One sees 
that for the hole-doped system, the tendency for holes to cluster together is very small - 
this is consistent with the binding energy study of Riera and Young []58| (if one uses their 
definition for the binding energy one finds results very similar to those of Table II). However, 
for the electron-doped material one finds that there is both a stronger binding of a single 
pair of holes, as well as a very low interaction energy, at least relative to the hole-doped 
material, as the number of carriers increases. This shows that the mobile electrons have a 
much greater propensity to cluster together than do the mobile holes. 

Summarizing the numerical work of the last two sections: (i) The single-hole quasipar- 
ticles of the hole-doped system are Shraiman-Siggia dipolar quasiparticles. They tend to 
remain far from each other, and are thus consistent with some form of superposition of these 
many-body quantum states in the low-doping limit |18| . (ii) The single-carrier ground state 
of the electron-doped material is different than that of the hole-doped system, viz. it is a 
(tt, 0) state as opposed to a (^, |) state, and a non-zero density of the former tend to exhibit 
a much stronger tendency to cluster together than do the holes of the hole-doped material. 
The spin distortions for these quasiparticles are long (short) ranged for the hole (electron) 
doped systems. 
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Considering the data contained in Fig. 6 and Table II, we have reason to believe that for 
the electron-doped materials the carriers have a strong tendency towards phase separation — 
certainly the tendency is much stronger than in the hole-doped system. However, since Table 
II shows that there are not purely attractive interactions between the mobile electron carriers, 
this tendency will probably lead to the presence of low-energy excited states corresponding 
to fluctuations into a locally phase separated state. Then, perhaps these phase separated 
states are the cause of the mid-gap states seen in the photoemission work [p6fl , as was 
suggested elsewhere p4| . 

As espoused by Emery and Kivelson []34}| , phase separation in the large hopping limit is 
likely to be frustrated, the frustration being due to the presence of long-ranged Coulomb 
interactions. We have investigated the effect of including a near-neighbour Coulomb re- 
pulsion between carriers on neighbouring sites, and find that no qualitative changes take 
place until very large Coulomb interactions (say of order much larger than t) are present. 
However, we do not know how long-ranged Coulomb interactions will effect our conclusions, 
and this is presently being studied. Further, Trugman noted that tendencies towards 
the binding of a pair of holes can be greatly reduced by next-nearest neighbour three-site 
spin-dependent hopping. While this form of hopping will be of much lower order than the 
next-nearest neighbour hopping that we are studying (since we begin with a three-band, 
not a one-band (as does Trugman), description), and thus we do not expect any qualitative 
changes in our analysis, we feel it appropriate to see in what quantitative ways our observed 
tendency towards a phase separated state in electron-doped planes changes with such a 
hopping term present. Thus, we are also presently studying the effects of such terms. 

Consequently, we propose that the commensurate antiferromagnetic spin background of 
the electron-doped materials do not have incommensurate correlations introduced into them 
because (i) the quasiparticles introduced into the weakly doped planes are not Shraiman- 
Siggia dipolar quasiparticles, and (ii) the mobile electrons found in these materials exhibit 
a form of phase separation, and no incommensurate correlations are expected for such a 
state of matter. We have completed a study of the electron/hole momentum distribution 
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functions, similar to those presented in Ref. ||18|| , including an evaluation of the density of 



states for both the hole and electron-doped compounds, and this work (to be presented 
later) further strengthens our assertion that some form of carrier clustering is taking place 
in the electron-doped materials. 

V. DISCUSSION: 

Our numerical results and semiclassical studies have led to the proposal that the single- 
hole ground states in the hole and electron-doped materials are different, and that this has 
a profound effect on (i) their magnetic phase diagrams, and (ii) the tendency towards phase 
separation that these carriers display. The Shraiman-Siggia dipolar quasiparticle picture 
has been shown p8| , p9| , |3T| to be consistent with many different experiments studying the 



hole-doped LaSrCuO system, and it will be interesting to see if our (jr, 0) quasiparticle 
picture of phase separation in the NdCeCuO can be used to explain more than just the 
ARPES results associated with the mid-gap states found in ARPES studies | f26f . 

One other experiment that may possibly find an explanation in the frustrated phase 
separation scenario for the electron-doped material involves the resistivity vs. temperature 
results of Hikada et al. fl33| . To be specific, if the mobile electrons tend to cluster together 
for sufficiently long periods of time, the dominant interaction between carriers is electronic 
as opposed to magnetic, and thus a fermi-liquid-like p ~ T 2 is not that surprising. Of 
course, such a qualitative explanation requires rigorous calculations to be carried out before 
it is to be taken seriously. 

We have addressed the magnetic and other material characteristics of the two 214 high 
T c systems. However, the approach of utilizing the t — t' — J system does not seem to be 
limited to these single-layer high T c materials. Recently, it has been proposed |59| that hole 



pockets appear in the YBaCuO material as one dopes away from the antiferromagnetic 
insulator, and that these pockets are centred at (f , f )• This is entirely consistent with the 
"small fermi surface" observed in ARPES studies of Y Ba^Cu^O^ z [BO]. In order to model 
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such a system one must include a t' parameter that rotates the fermi surface by 45° with 
respect to that predicted by the t — J model, and to this end we have used t — 1, t' ~ —.5, 
and J = .4. Then, we find (i) that indeed the single-hole ground state occurs at k = (n, n) 
(this was also found in the study of Ref. ||44||), and (ii) that pairs of such holes tend to 
stay as far apart as possible, with a hole-hole correlation function similar to that shown for 
the hole-doped material in Fig. 6. Thus, weakly interacting quasiparticles forming small 
hole pockets at (n, it) is consistent with our numerics, and are not expected to lead to any 
incommensurability ||59|| . Noting that both our semiclassical work (similar to that shown in 
Fig. 5) and the spin-wave analysis of Reiter predict that the spin distortions produced by 
(7r, 7r) quasiparticles are long ranged, it is very interesting that in the two systems in which 
the holes tend to stay as far apart as possible, the single-carrier quasiparticles produce 
long-ranged spin distortions, and we are presently exploring this coincidence. 
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FIGURES 

FIG. 1. A schematic of the phase diagrams, as a function of temperature and doping concen- 
trations away from half filling, for both the hole and electron-doped materials (after Ref. [36], 
fromRef. [4]). 

FIG. 2. The magnetic structure factor for the t — t' — J model given in Eq. (2.4) for 0, 1, 
2, 3 and 4 carriers; here these carriers are mobile holes. For this figure, t = +1, t' = —.3, and 
J = .4, which are appropriate for hole-doped LaSrCuO. The reciprocal lattice points are as 
follows: T = (0,0), X = (vr,0), and M=(vr,7r). 

FIG. 3. The magnetic structure factor for the t — t' — J model given in Eq. (2.4) for 0, 1, 2, 
3 and 4 carriers; here these carriers are mobile electrons. For this figure, t = —1, t' = +.3, and 
J = .4, which are appropriate for electron-doped NdCeCuO. 

FIG. 4. A picture of the equivalent sites that we allow to be distorted away from their ordered 
Neel configuration in our calculation of the semiclassical wave function for the t — t' — J model for 
one carrier. The clusters used in such calculations are infinite in extent. 

FIG. 5. The energies of the k = (ir,0) ground state for the electron-doped system 
(t = —1, t' = +.3, J = .4) and the k = (§,§) ground state for the hole-doped sys- 
tem (t = +1, t' = —.3, J = .4) obtained via application of the variational principle to a 
Shraiman-Siggia semiclassical wave function. These energies are found by allowing all spins up to 
distance d, given in units of the lattice constant a, to be distorted away from their antiferromag- 
netically aligned state. 

FIG. 6. The carrier-carrier correlation function, defined and described in Eq. (6.1) of Ref. [18], 
for the hole (t' = —.3) and electron (t' = +.3) doped systems, evaluated for the doubly doped 
state. It is seen that the holes tend to stay as far apart as possible, whereas the mobile electrons 
tend to cluster together at short distances. 
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TABLES 

TABLE I. Single carrier minimum energies for both hole and electron singly-doped planes, 
described by Eq. (2.4), for the allowed wave vectors of our non-square 24-site cluster. All energies 
are in units of t. 



k 


Energy 


Energy 




(Hole doped: t = 1, t' = -.3) 


(Electron doped: t = — 1, t' = +.3) 


(0,0) 


-12.824 


-12.153 


/ 7T 7T\ 
1/11/1/ 


-12.961 


-12.338 


/ 7T 7T \ 
v A 1 A 7 


-13.207 


-12.686 


/ 37T 37T \ 

\ A 1 A ) 

\ 4 / 4 / 


-13.063 


-12.375 


(7T,7r) 


-12.853 


-12.177 


/ — TT TT \ 

\ 3 ' 3 > 


-12.997 


-12.471 


I -2tt 2tt \ 
V 3 ' 3 > 


-13.127 


-12.363 


1 — TT 7lT \ 

V 12 ' 12/ 


-12.896 


-12.746 


/ — 5n Hit \ 
V 12 ' 12 1 


-12.980 


-12.679 


ITT 5TT \ 

V6' 6 > 


-12.821 


-13.037 
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TABLE II. Interaction energies, denned in Eq. (4.1), for 2, 3 and 4 carriers in the t — t' — J 
model, using representative material parameters for the hole and electron-doped compounds based 
on Eq. (2.4). As long as the single-hole ground states are those that we predict for the hole and 
electron-doped systems, these interaction energies do not sensitively depend on the ratio of | t'/t \ . 





Hole Doped 


Electron Doped 




t = 1, t' = -.3 


t = -1, t' = +.3 


E 2 / 


-.12 


-.21 


E37 


+.45 


+.12 


E4I 


+1.5 


+.26 
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